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We consider the electron-positron plasma generation processes in the magnetospheres 
of magnetars — neutron stars with strong surface magnetic fields, Bq ^ 10^'^— 10^^ G. We 
show that the photon splitting in a magnetic field, which is effective at large field strengths, 
does not lead to the suppression of plasma multiplication, but manifests itself in a high 
polarization of 7-ray photons. A high magnetic field strength does not give rise to the sec- 
ond generation of particles produced by synchrotron photons. However, the density of the 
first-generation particles produced by curvature photons in the magnetospheres of magne- 
tars can exceed the density of the same particles in the magnetospheres of ordinary radio 
pulsars. The plasma generation inefficiency can be attributed only to slow magnetar rota- 
tion, which causes the energy range of the produced particles to narrow. We have found a 
boundary in the P—P diagram that defines the plasma generation threshold in a magnetar 
magnetosphere . 
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INTRODUCTION 

Magnetars are neutron stars with superstrong surface magnetic fields, 
Bq ~ 10^^—10^^ G. This is two or three orders of magnitude higher than the field 
strength typical of radio pulsars, 10^^ G. Such magnetic fields for magnetars follow from 
the standard estimate of Bq ~ 2(PA-i5)^/^ 10^^ G that is used for radio pulsars and that 
is based on the fact that the neutron star rotation energy is spent on the generation of 
magnetodipole radiation, which owes its origin to the magnetic field rotation. Just as radio 
pulsars, magnetars spin down, P > 0, but with a much higher rate, P ~ 10^^°— 10^^^, 
while the typical spindown rate for radio pulsars is P ^ 10~^^,P_i5 — 1- The rotation 
periods P of magnetars are also longer than those of radio pulsars and lie within the range 
P ~ 5—10 s. The application of the magnetodipole formula is more justified for magnetars 
than for radio pulsars, whose magnetospheres must be filled with a dense electron-positron 
plasma that shields the magnetodipole radiation (Beskin et al. 1993). In addition, the 
observation of an absorption line in the X-ray spectrum of the magnetar SGR 1806-20 
(Ibrahim et al. 2002), which is interpreted as a proton cyclotron resonance near the stellar 
surface, yields the same estimate of the magnetic field strength for magnetars. 

The term "magnetar" was introduced by Duncan and Thompson (1992) for soft gamma 
repeaters (SGRs). At present, two groups of objects are classified as magnetars: SGRs and 
anomalous X-ray pulsars (AXPs), which currently number 5 and 9, respectively. SGRs yield 
gamma-ray bursts with energies up to 3.5 x 10^^ erg, which was observed on December 27, 
2004, in SGR 1806-20; AXPs also exhibit bursts but in the X-ray energy range. SGRs and 
AXPs have similar periods and magnetic fields, but P for SGRs is slightly larger, P ~ 10"^". 
Therefore, SGRs may be considered to be younger objects, with ages of ~ 10^ yr, than 
AXPs, whose ages are ~ 10^—10^ years. 

It is important to note that the magnetic fields of magnetars also determine their energy 
activity. Indeed, the energy stored in the magnetic field of a magnetar, 

exceeds the neutron star rotation energy. 

This is not the case for radio pulsars: Wr ^ Wb- In addition, the X-ray luminos- 
ity of AXPs, Lx — 10^^—10'^® erg s~^, is much higher than the rotation energy being lost, 
Efi = AiT^IP^^P ~ 10'^^— 10^^ erg s~^. This means that the energy being released is not the 
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neutron star rotation energy, as is the case for radio pulsars, but the magnetic field energy 
stored in the star. How the magnetic field energy of magnetars is released and where such 
strong fields in them originate from still remains a puzzle. In addition, it is unclear why 
the spindown of the star is still observed if its rotation is not the energy source. Strong 
magnetic fields can be generated by magnetic dynamo mechanisms (Duncan and Thompson 
1992; Bonanno et al. 2006) and can result from the penetration of the radiation field into 
the neutron star crust with a high but finite conductivity (Istomin 2005). 

The absence of radio emission from magnetars in the standard range of frequencies at 
which radio pulsars are observed, v ~ 0.5—3 GHz, suggests that there is no dense plasma 
in the magnetospheres of magnetars whose fiows would generate radio emission, as is the 
case in the magnetospheres of radio pulsars. However, low-frequency (60—110 MHz) radio 
observations at the Pushchino Radio Astronomy Observatory revealed emission from some 
magnetars (Malofeev et al. 2004, 2005), although its characteristics differ from those of radio 
pulsars. 

The goal of this paper is to investigate the generation of electrons and positrons in 
the magnetospheres of neutron stars with superstrong magnetic fields which are present in 
magnetars. We wish to ascertain why no dense plasma is generated in magnetars and what 
suppresses its effective generation. 

For magnetars, the following critical magnetic field strength is a characteristic one: 

2 S 

B, = — ^ ^ 4.414 X 10^^ G, (1) 
eh 

where m is the electron mass, e is the electron charge, c is the speed of light, and h is the 
Planck constant. When the magnetic field strength becomes comparable to (1), the quantum 
nature of physical processes begins to manifest itself. If we consider the motion of an electron 
in such a field, then only its momentum component parallel to the magnetic field direction 
remains continuous. The transverse momentum component is quantized in such a way that 
the expression for the total particle energy is 

el = l+p^ + 2nB, (2) 

where p is the longitudinal particle momentum component and n = 0, 1, ... is the number 
of the Landau level occupied by the particle. In what follows, we will measure the magnetic 
field strength in units of (1) and take the electron mass m, the Compton wavelength of the 
electron X = h/mc ~ 3.86 x 10^^^ cm, and its ratio to the speed of light A/c ~ 1.29 x 10~^^ s 
as the units of mass, length, and time, respectively. Formally, this means that h = X = c = 1. 
Thus, the length, velocity, time, energy, and wave vector will be measured in units of X, c, 
A/c, mc^, and 1/A, respectively. 
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Two distinctly different regions are known to be identified in the magnetospheres of 
rotating neutron stars: tlie regions of open and closed magnetic field lines. In the region of 
closed field lines, the particles rotate synchronously with the field because of high plasma 
conductivity and magnetic field freezing-in. However, along the open field lines, the particles 
can move freely and escape from the neutron star magnetosphere. Continuous plasma outflow 
from the magnetosphere requires the presence of electron-positron pair generation processes 
compensating for it. 

In this paper, we consider the generation of relativistic electrons and positrons by pho- 
tons in the superstrong magnetic fields of magnetars, the photon splitting in such fields, 
and its influence on the pair generation. By taking this into account, we derive expressions 
for the particle distribution function and density and find a criterion for effective plasma 
generation. 



KINETIC EQUATION FOR RELATIVISTIC PARTICLES 

The electron-positron plasma generation processes in the magnetosphere of a neutron 
star are described by the kinetic equation 

Here, F is the particle distribution function in longitudinal momentum p or in longitudinal 
Lorentz factor 7, which in dimensionless units coincides with the energy (2) at n = 0, 
V is the particle velocity, b = B/i? is a unit vector along the magnetic field, Q{F) is the 
operator that describes the generation of particles by high-energy photons, and S{F) is the 
operator that describes the scattering of particles through the emission of photons by them. 

A particle with a Lorentz factor 7 moving along a field line with a radius of curvature p 
at a given point emits photons called curvature ones. The emission probability P^\''^,k) 
of a curvature photon with energy k per unit time is given by the expression (Sokolov and 
Ternov 1983) 

p'i'M-^^^.i'^l (4) 



2^3 717^ 

where a = /he = 1/137 is the fine-structure constant, P^\'~f, k) is the emission probability 
of a 1 1 -polarized photon whose electric field vector lies in the plane passing through the 
vectors k and B, and P^\'^,k) is the emission probability of a ±-polarized photon whose 
electric field vector is orthogonal to this plane. The characteristic energy k^ of the emitted 
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curvature photons is 

The functions <f±{x) are defined as follows: 

oo 

'f±{x) = J K^i?Xy) dy ± K2/3{x), 

X 

where K^{x) is the Macdonald function. Using asymptotic representations of the Macdonald 
function, we can easily derive asymptotic expressions for the functions and ^P-{x) that 

are useful in the subsequent analysis: 

2-1/33 r[2\„2/3 fora;-^0 

Mx)={ (6) 
I — e ^ for a; — >• 00; 

2-i/3r('-^ a;-2/3 forx^O 

^-i^)={^^^^^ (7) 

X ' e tor X 00. 

12 



Note that Eq. (4) is valid if the classical theory of synchrotron radiation is applicable. 
The criterion for its applicability is 

— < 1. (8) 
P 

Since the Lorentz factor of the particles 7 and the radius of curvature of the field lines p in 
the magnetospheres of neutron stars are typically 7 < 10^ and p > 10^^, in dimensionless 
units (note that in these units, 1000 km ^ 0.26 x 10^^), condition (8) is always satisfied. 
The satisfaction of criterion (8) means that the characteristic energy of the emitted photons 
accounts for only a small fraction of the initial particle energy (see (5)). 

Klepikov (1954) was the first to show that a single-photon generation of elec- 
tron-positron pairs is possible in a magnetic field. Let a photon with energy k move at 
an angle x the direction of magnetic field B. The photon absorption coefficient due to 
electron-positron pair creation in the case of a weak magnetic field, -B ^ 1, for two directions 
of photon polarization is then (Tsai and Erber 1974) 

/ 8 1 \ 

fi\\,±iB,k,x) = bi±aB sinxexp ( - --^z—. ), (9) 

V S Bk smy/ 
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where 

\/3 \/3 
6ll = ^^0.31, 6x = ^~0.15. 

Equation (9) is applicable when Bk sinx <^ 1- We should also take into account the fact 
that the pair creation is possible only when a photon with energy k propagating along the 
tangent to the magnetic field line at the point of emission reaches some threshold angle Xt 
at which the transverse component of the photon wave vector is equal to 2, 

ksmxt = 2. (10) 

If, alternatively, we consider the case of a strong magnetic field, S > 1, characteristic 
of magnetars, then we should take into account the fact that the electron-positron pairs are 
generated in such a way that either both particles of the produced pair are at the zeroth 
Landau level if the pair-generating photon is || -polarized or one particle of the pair is at the 
zeroth Landau level and the other particle is at the first Landau level if the pair-generating 
photon is _L-polarized (Semionova and Leahy 2001). In the former case, the photon absorp- 
tion coefficient is 

2aB fk"^ sin^x \^^^'^ / Psiri^x 



and the threshold angle Xt meets to the same condition (10). In the latter case, the absorption 
coefficient k and the angle Xt are 

K^{B, k, x) = —r^. 1 + pIi + SoBi) exp — ^ , 12 

Poi^ smx ^ ' V 2B / 



where 



A; sin Xi = 1 + Vl + 25, 

= ~A 1 - 5 + . 2 

\ 4 k'^ sm X 
^0 = ^ -B + . o , £1 = ^ + B + 



k"^ sin x) \ 4 k"^ sin x 

Let us now find the form of operator Q in the kinetic equation (3), which describes the 
generation of electron-positron pairs. Let qo{k) be the number of photons emitted per unit 
volume per unit time. The opening angle of the emission cone is inversely proportional to the 
particle Lorentz factor 7. Since the energy of the particles is fairly high, all of the emission 
may be assumed to be directed along their motion. The distribution function fo{k) of the 
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emitted photons will then be proportional to the radius of curvature po of the magnetic field 
line at the point of emission: 

m) = poqoik). (13) 

The function fo{k) is normalized in such a way that fo{k) dk gives the number of photons 
with an energy in interval dk about k. Obviously, if there is no photon absorption, then 
the distribution function is conserved along the direction of emission. If, alternatively, the 
electron-positron pair generation takes place, then the photon distribution function f{k) 
at the point under consideration is smaller and can be expressed in terms of the photon 
distribution function fo{k) at the point of emission and the absorption coefficient k as 

fik) = foik)exp(^-J K{B{l'),k,xil'))dl'), (14) 
h 

where / is the distance from the point under consideration and the point of emission and It is 
the so-called threshold distance at which the electron-positron pair generation begins. 

Assuming that the magnetic field in the neutron star magnetosphere in the polar region 
of interest near the surface is a dipole one, let us introduce the cylindrical (r, z) coordinates 
in such a way that r denotes the distance from the dipole axis to the point of observation 
and z denotes the distance from the dipole center coincident with the stellar center to a 
given point along the dipole axis. Expanding the well-known expression for the shape of a 
dipole magnetic field line = Rsiia^ 9p specified in polar {rp,9p) coordinates at small polar 
angles 9p, we find that 

' (15) 



where R> Rc is the distance from the stellar center at which a given field line will cross the 
z = plane when extended as a dipole one, Rc = c/Q is the light-cylinder radius, and Q is 
the pulsar rotation rate. We derive an expression for the radius of curvature of the field line 
from Eq. (15): 

p = ^V^. (16) 

Assuming that z = Rg and R = Rc, where Rg ~ 10 km is the neutron star radius and 
i?c ~ 3 X 10^ km is the characteristic light-cylinder radius for a pulsar with a period P = 6 s, 
we obtain the radius of curvature p > 10^ km mentioned above. Note that Eqs. (15) and 
(16) are valid to terms 0{z/Rc), which may be disregarded with a sufficient accuracy in the 
conditions under consideration. 



Using Eqs. (15) and (16), we can easily obtain the dependence of the angle x between 
the directions of magnetic field B and photon wave vector k on the distance / traversed by 
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the photon: 

x{l) = -ttV' (1^) 

Po 1 + l/Zo 

where po is the radius of curvature of the magnetic field hne at the point of emission and Zq 
is the distance to this point from the dipole center. Since we used the fact that the angle 
between the directions of the magnetic field and the photon wave vector is small, we can 
replace sinx by x- Below, we will assume that / <^ zq. This means that 

X(/)«Xmax=-, (18) 

Po 

where Xmax is the maximum possible angle between the photon wave vector and the magnetic 
field direction. When criterion (18) is met, we may ignore the dependence of the magnetic 
field strength B on the distance traversed by the photon when integrating in Eq. (14) and 
assume (see (17)) that 

X{1) = -. (19) 
Po 

The number of electron-positron pairs produced per unit time per unit volume at the 
point of absorption is 

oo 

dN 



dt 



fik)K{B,k,x)dk, (20) 



kt 



where, as in the case of Eq. (10), the threshold energy of a photon capable of producing 
electron-positron pairs, kt = 2/x, is defined. 

According to Beskin (1982), the photon energy for a weak field in the case of elec- 
tron-positron pair generation is distributed almost evenly between the electron and the 
positron. It thus follows that the energy of the produced particle is uniquely related to the 
angle between the photon wave vector and the magnetic field direction: 

7 = -, (21) 
X 

where 7 is the Lorentz factor of the particle after it passes to the zeroth Landau level. For 
a strong field, S ^ 1, as we will see below, the electron-positron pair generation threshold 
matches (10), since all photons will be ||-polarized due to the photon splitting in a strong 
magnetic field when this threshold is reached. In addition, the pair generation in this case 
will take place at the threshold point at angles (x — Xt)/Xt ^ 1- Therefore, the magnitude 
of the pair particle momentum in the reference frame where k ± B is close to zero. Hence, 
using the inverse Lorenz transformation, we can easily find that equality (21) is also valid 
for > 1. 
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Thus, we find from (13), (14), (20), and (21) tliat tlie operator Q describing tlie particle 
source in tlie kinetic equation (3) is 

oo X 

Q = Poj ^^^{^^~) J dkqoik) K{B,k,x)exp{-po j K{B,k,x') dx'), (22) 

2/x 2/k 

where B is the magnetic field strength at the point of absorption and the absorption coeffi- 
cient K is defined by Eq. (11) for the strong magnetic fields of magnetars or by Eq. (9) for 
the weak magnetic fields of ordinary pulsars. In Eq. (14), we passed from integration over 
the distance I' to integration over the angle x' using (19). 



Gurevich and Istomin (1985) gave the form of the operator 5* describing the scattering 
of particles through the emission of curvature photons by them: 

S=— — 
3p^ ^7 



\ 55 di^'F) 

7 -T + 



(23) 



32^3 p <97 

Equation (23) is valid when criterion (8) is met. These authors pointed out that the oper- 
ator 5* gives only a small contribution compared to the operator Q. Therefore, it may be 
disregarded when the electron-positron pair generation processes are considered. 

The last term on the left-hand side of Eq. (3) is responsible for the particle acceleration. 
We will assume that this acceleration takes place in a small region above the neutron star 
surface where the so-called inner gap is formed (Ruderman and Sutherland 1975). A signif- 
icant accelerating electric field exists in this region, but outside it this field is shielded by 
the generated plasma and no particle acceleration takes place. It thus follows that we may 
disregard the term under consideration and take into account the particle acceleration as a 
boundary condition at the neutron star surface for the particle distribution function, 

F(7) = iVo5(7-7o), (24) 

where Nq is the electron-positron plasma density at the neutron star surface and 70 — 10"^ 
is the Lorentz factor that the particles acquire as a result of their acceleration in the inner 
gap. 

As regards the outer gap (Cheng and Ruderman 1977), if it emerges in the magneto- 
sphere of a magnetar, then the electron-positron plasma generation processes in it do not 
differ from those in an ordinary radio pulsar, since it is located far from the surface where 
the magnetic field is weak 

Thus, the kinetic equation (3) in the case under consideration transforms to 

^ + div(Fv) = Q, (25) 

where the operator Q is defined by Eq. (22). 



PHOTON SPLITTING IN A STRONG MAGNETIC FIELD 

The photon sphtting processes become important in the strong magnetic field of a mag- 
netar. It should be noted that the splitting probabilities are very sensitive to the photon 
polarization. Based on the CP invariance and the energy-momentum conservation law with 
vacuum polarization, Adler (1971) showed that -L— ^|||| decay is crucial in moderately strong 
fields, -B ^ 1, and at low photon energies. Before the appearance of the paper by Usov 
(2002), the case of ||^-L|| and ±— decays in the calculations of magnetar magneto- 
spheres, in general, was not ruled out (Harding et al. 1997; Baring and Harding 2001). How- 
ever, Usov (2002) showed that Adler's selection rules are valid in fields of arbitrary strengths. 
Thus, the decay of || -polarized photons below the electron-positron pair generation thresh- 
old (10) is strictly forbidden. Hence, it would be incorrect to consider a polarization-averaged 
photon distribution function and to use a polarization-averaged attenuation coefficient in a 
strong field. 

Let /|| (k) and f±{k) be the photon distribution functions in energy k for two directions of 
polarization. We assume that the splitting of only _L-polarized photons is possible. However, 
the photon need not decay into two photons with equal energies. We will describe the 
splitting probability of a ±-polarized photon with energy ki + k2 into two ||-polarized photons 
with energies ki and k2 per unit time by a function w{ki,k2)- 

The system of equations that describes the evolution of the distribution functions for 
photons of two polarizations is 

k 

f±{k) = -\fdk) j w{k-k\k')dk', 



oo 

f\\{k) = j w{k,k')f^{k + k')dk'. 



The first equation in (26) describes the splitting of a ±-polarized photon with energy k 
into two ||-polarized photons with energies k' and k — k'. The second equation describes 
the generation of a || -polarized photon with energy k through the decay of a photon with 
a higher energy k + k' and the detachment of a photon with energy k'. We immediately 
note that time differentiation for the photon distribution functions is equivalent to distance 
differentiation, since we assume that c = 1. In this case, the splitting probabilities per unit 
time and per unit length of the distance traversed by the photon are equal. 

It follows from the first equation in (26) that the total photon splitting probability per 



(26) 
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unit time is 

k 

Wik) = ^ J w{k- k', k') dk'. (27) 



Let a photon with a wave vector k be emitted by relativistic particles propagating along 
the magnetic field lines at some point in the direction of their motion. Baring and Harding 
(1997) provide an approximation of the total photon splitting probability per unit length in 
the limit i? ^ 1 as a function of the photon energy k with an accuracy of ~ 1.5%: 

W{kX) = ^?f^ + — C - — C + — exp(5C)l . (28) 

^ '^^ TT^ k V271 124^ 247^ 297^ 3461 ^7 ^ ' 

Here, we introduced a variable = kx/^, where x is the angle between the vectors k and 
B related to the distance traversed by the photon by Eq. (19). The equation that describes 
the change in the distribution function of _L-polarized photons is 

C 

/±(fc,C) = /±(A:,0)exp(-^ I W{k,C)dC), (29) 



where f±{k, 0) is the distribution function of ±-polarized photons at the point of emission. As 
the photon moves, the variable ( increases and becomes equal to unity when the threshold of 
electron-positron pair generation by || -polarized photons (10) is reached. Directly integrating 
Eq. (28) and substituting ( = I yields 

/±(A:,l) = /±(A;,0)exp(-3.57x 10-3^0). (30) 

Thus, we may assume that all _L-polarized photons had decayed before the threshold 
of pair generation by || -polarized photons was reached, when the the expression in the ex- 
ponent > 1. Therefore, the inequality 0.56 x 10~^ po/k"^ > 1 must hold. Substituting the 
characteristic radius of curvature of the magnetic field lines, p > 10^^, into this inequality 
yields 

k < 10'. (31) 

Note that inequality (31) always holds for typical conditions in the magnetosphere 
of a magnetar, since the maximum energy of the curvature photons has a characteristic 
value of A; ~ 10^ for a particle Lorentz factor 7 ~ 10^ (see (5)). Above this energy, the 
curvature-photon distribution function decreases exponentially (see asymptotics (6) and (7)). 

We see that as a result of the splitting of ±-polarized photons for magnetars, all photons 
with be 1 1 -polarized when the threshold angle (10) is reached. 
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Baier et al. (1996) showed that the total photon sphtting probabihty increases as the 
electron-positron generation threshold C = 1 is approached. In this case, the differential 
photon splitting probability near the threshold is almost constant, i.e., the photon can decay 
into any two lower-energy photons almost equiprobably, provided that the sum of their 
energies is equal to the energy of the original photon. If, alternatively, ( 1, then, in 
general, this is not true. However, let us turn our attention to the approximation of the total 
splitting probability (28). The photon splitting probability increases with (; at small (, the 
probability W{k, () oc It thus follows that the splitting will take place mainly near the 
pair generation threshold C = 1- At small (, the total splitting probability is so low that 
the specific form of the differential probability w{ki,k2) is unimportant, provided that its 
integration according to (27) yields an exact value for the total probability W{k). Therefore, 
we will assume below that the photon decays equiprobably into two lower-energy photons 
even at small (. 

If the _L-polarized photons split equiprobably into any two || -polarized photons, then 

W{ki,k2) =w{ki + k2). (32) 
Given Eq. (32), the system of equations (26) transforms to 

f±{k) = -'^w{k)h{k), 

7 (33) 
/|,(fc) = J w{k')f^{k')dk'. 

k 

Expressing the product w{k)f±{k) from the first equation in (33) and substituting it into 
the second equation yields 

oo 

^/ii(fc,C) .5 r dk' 

k 

Here, we used the fact that d/dt = d/dl = [k/{2po)] d/d(. We are interested in the value 
of the distribution function for ||-polarized photons a.t ( = 1. Below, we will see that the 
generation of electron-positron pairs by curvature photons in the strong magnetic field of 
a magnetar actually takes place at ( = 1. It should be noted that the splitting at ( = 1 
has already finished and there are virtually no _L-polarized photons (see Eqs. (30) and (31)). 
Therefore, we should set f±{k,l) = 0. Integrating (34) by taking this circumstance into 
account, we find that 

oo 

/dk' 
h{k',0) — . (35) 
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We see from Eq. (35) that for equiprobable splitting of ±-polarized photons (see (32)), the 
asymptotic expression for the distribution function of || -polarized photons does not depend 
on the specific form of w{k). 



THE ELECTRON-POSITRON PLASMA GENERATION 

Let us derive the specific form of the operator Q that describes the generation of rel- 
ativistic electrons and positrons by high-energy photons for a weak magnetic field, -B ^ 1, 
characteristic of ordinary pulsars, and for a strong magnetic field, B > 1, characteristic of 
magnet ars. 

Consider the case of i? ^ 1. If we pass from the variable x' to /i' = 2/ {kx') and from 
the variable k to fi = 2/{kx) in Eq. (22) for the operator Q, then it will transform to 

1 1 
Q = 2ahBpoj dxS[^--) I ^go(— )exp(-a/z-a6poxV/ ^SV"'^'), (36) 

A* 

where by the coefficient b we will mean either b\\ or b± (see (9)). When criterion (18) is met, 
we may assume that B' does not depend on the variable /i' and coincides in value with the 
quantity B. Integrating the expression in the exponent by parts and retaining the principal 
term (which is legitimate, because, as we will see below, the integrand differs noticeably from 
zero only at ^ 1), we obtain 

a 

<5 = / ^ go (27^) exp (^-y - . (37) 







Here, we use the notation 



abpoB 4 
P=— ' « = ^- (38) 



Let us determine the maximum of the exponent in Eq. (37). Assuming that this maxi- 
mum is reached aiy = A ^ 1, we find that A must be a solution of the equation A = Aq— In A, 
where Aq = In p. At large A, the solution of the equation is given with a sufficient accuracy 
by 

A^Ao-lnAo. (39) 

For typical radius of curvature of the field lines po ^ 10^^, magnetic field B ~ 0.01—0.1, 
and energy of the produced particles 7 < 10"^, A ~ 15—20. Since the coefficient b appears 
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under the logarithm, it makes no sense to distinguish Ay and Aj^ for different directions of 
photon polarization. Below, we will assume that b = + b±^)/2 ^ 0.23. 

Expanding the exponent at y = A gives exp(— y — pe^^/y) ~ exp(— A — {y — A)^/2) 
at A ^ 1. We see that this function has a sharp peak aX y = A and that its variance o"^ = 1. 
Therefore, the peak width is equal to unity and does not depend on A. 

As we see from (37), if a < A, then Q = 0. If, alternatively, a > A, then for a weak 
magnetic field, 

Q = ^go(27j), B«l. (40) 
We used the fact that at A 1, the integral 

oo 

f dy e-y ( e-y\ 1 

J y y ^ 2/ ^ A 



where A is defined by (39). 

Let us now consider the case of i? > 1. Using (22) and (11), we can easily obtain 

7r/2 

^ ^ apo5 f ^^cosr^exp^- ^ )qo(^^) 



where cost] = 2/{kx) and the angle x is related to 7 by Eq. (21). Since the condition 
^ 1 is always satisfied, the integrand differs noticeably from zero only at the lower 
limit of integration. Expanding the integrand at 77 = and performing integration from zero 
to infinity, we find that for a strong magnetic field, 

Q = 2go(27), B>1. (41) 

The physical meaning of Eqs. (40) and (41) is simple. At 5 ^ 1, an emitted photon with 
energy k will traverse the distance If = It a/ A, where It = 2po/^, and will then immediately 
generate an electron-positron pair at this point. Thus, // is the photon mean free path in a 
magnetic field. If i? > 1, then If = It, i.e., for the strong magnetic fields of magnetars, the 
generation of pairs by photons takes place immediately after the threshold has been reached. 

When condition (18) is satisfied, the value of qo{k) at the point of photon emission is 
equal to the value of q{k) at the point of photon absorption, i.e., qo{k) ^ q{k). Let us derive 
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an expression for the function q{k). For a weak magnetic field, 

q{k) = q^'\k) + q^'\k), (42) 

where the functions q^^\k) and q^'^\k) give the numbers of curvature photons and syn- 
chrophotons emitted by a unit volume per unit time. We will distinguish the functions 
q^\k) for ||-polarized photons and q^l\k) for ±-polarized photons: 

oo 

qt\k) = I Pf^(7,fc)Fo(7)rf7. (43) 



The lower limit in Eq. (43) may be assumed to be zero, because criterion (8) is met. 

To calculate the specific form of the function q^'^^k) for i? ^ 1 and i? > 1, we must 
know the form of the distribution function for the primary particles -^0(7) after their accel- 
eration in the inner gap. When the boundary condition (24) is satisfied, the primary particle 
distribution function is 

^o(7) = iVo|^5(7-7o), (44) 

-Do 

where B is the magnetic field strength at the point under consideration, Bq and No are the 
magnetic field strength and the plasma density at the stellar surface. We will assume that Nq 
is proportional to the Goldreich- Julian density |A''^'^| = QBQ\cos9\/{2'7ice) at the stellar 
surface, where 6 is the inclination of the magnetic dipole axis to the neutron star rotation 
axis, with the proportionality factor being determined by the total longitudinal current j\\ 
flowing in the magnetosphere, so that A^o = ^ol-^f ''I; where io = 27rj\\/[QBo \ cos6\) ^ 1. 

For the weak magnetic fields of ordinary pulsars, 

,«(l.)=,?W + ^'W = ^4#*'(f^). (45) 

V3 TT 7o ^0 ^ 37^; / 

where the function ip{x) is defined as 

00 

fix) = ^ = / K5/3{y)dy (46) 

X 

and has the asymptotics 




- 16 - 



If, alternatively, B > 1, then the photon decay processes are important. Therefore, Eq. (35) 
must be used to find the form of the function q^'^\k). Then, 

<l''\k) = ^-,'^9m). B>1, (48) 
where the function g{x) is defined as 

oo 

g{x) = M^) + 2 [ v4y)-- (49) 

J y 



Its asymptotics are 



22/33 r( x-2/3, X^O 

9i^) = < p^^^^ (50) 




X 



For i? ^ 1, no synchrophotons are emitted, since the particles are produced straight at 
the zeroth Landau level. Therefore, we should set q^^\k) = in Eq. (42). We see that only the 
first generation of particles produced by curvature photons exists in the magnetosphere of a 
magnetar. If, alternatively, <C 1, as in the magnetosphere of a pulsar, then the contribution 
from synchrophotons is significant and q^^\k) is nonzero; therefore, the second generation of 
particles can be produced by synchrophotons. For comparison with the magnetosphere of a 
magnetar, we will need only the first-generation particle distribution function. 

Given the conditions 7^1 and divB = 0, the kinetic equation (25) for the steady-state 
case transforms to 

bVF{ = Q', (51) 

where F[ = Fi/ B and Q' = Q/B are, respectively, the distribution function for the first gen- 
eration of particles and the particle source divided by the magnetic field strength. If i? ^ 1, 
then Q can be determined from (40) and (45). So, Eq. (51) takes the form 

d 3r d \ ^,(^) 2 a Nq a (^"^p a 



We see from Eq. (47) for the asymptotics of v^(x) that the first-generation particle 
distribution function F[{'-)) for a weak magnetic field differs noticeably from zero only 
at 7 < 7^x = 37QA/(4ap). At 7 > 7max5 the function -^{(7) decreases exponentially and 
gives no contribution to the total plasma density. We will be interested in the range 7 < 7^L 
In this case, the solution to Eq. (52) is 

5/3 



- 1 



(53) 
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5 « 1, 7 < 7^1 = 37o'A/(4ap). 

Here, Rg is the neutron star radius, Oq = 4/(3-Bo), -Bo is the surface magnetic field, 
and R = z^/r"^ is the integral of motion that is conserved along a magnetic field line. We set 
the boundary condition -^{(7) = at z = Rs, because the curvature photons producing the 
first generation of particles just begin to appear near the stellar surface. 

Let us now consider the case of a strong magnetic field, B > 1. Using Eqs. (41) and (48), 
we obtain 

d 3r d\ ^,(,) _ 1 a A^o „f^1P\ 



B > 1. 



(54) 



Just as in the previous case, we will be interested in the form of the distribution func- 
tion -F{(7) only at 7 < 7max = 37q/(4p) (see asymptotics (50)). Under this condition, 
Eq. (54) has the following solution: 



"StT 



r(-la^i?-l/3i?2/3^-2/3 



Bn 



2/3 



(55) 



B>1, 



7 < 7il = 37oV(4p). 



We see that the first-generation particle distribution functions for both weak and strong 
magnetic fields are F[('~f) cx 7"^/^ and have the same order of magnitude. Near the surface, 
the functions increase linearly in [z — Rg), with their ratio being 



F{(-)(7) _2/aoy/3^^ 
F/(^)r^^ 3VA/ 



The density ratio will be slightly different, because 7max depends on the magnetic field 
strength, and is smaller in a weak field than in a strong one. Integrating Eqs. (53) and (55) 
over the particle Lorentz factor from 7inin to 7inax yields 



' lOvr [3) Bo\a) ' ^'"^^ 



1 - 



7min 

Iw}' 
7max 



1/3- 



Z 

R. 



5/3 



- 1 



, (56) 



5 < 1, 



7min < 7 ^ 7maxi 



8ir Vi B„ • 



7min 
7max 



1/3- 



Z 

Rf. 



2/3 



- 1 



(57) 



7mm < 7 < 7 



is) 
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The minimum particle Lorentz factor ^ram is determined by the maximum angle Xmax 
(see (18)) and is equal to 7min = 1/Xmax — 100. Since the true values of the distribution 
function -F(7) = BF'{'~f) and the plasma density = BN', the ratios of the distribution 
functions and densities near the surface for strong and weak magnetic fields are 



iw) 



Oo 

A 

2 R 



1/3 R 



iw) 



5, 



(w) 



(^) 



(5^ 



(59) 



while 5^"^ and B'^"' are the surface magnetic field strengths for an 



where = 4/(35^"'^), while B^""^ and B^^^ 
ordinary pulsar and a magnetar, respectively. In calculating ratios (58) and (59), we used 
the conditions 7min -C 7max and {z — Rs)/Rs ^ 1. 

A factor that increases when 2; — > 00 appears in Eqs. (53), (55), (56), and (57). This is 
because we assumed condition (18) to be satisfied, i.e., we assumed that the mean free path of 
a curvature photon is small; therefore, the distribution functions and magnetic field strengths 
at the points of photon emission and absorption coincide. However, if x ~ Xmax = zq/pq, 
then the mean free path is large and the electron-positron pair generation takes place at 
distances much larger than zq. Therefore, the rate of particle generation decreases with 
increasing z and the distribution function -^{(7) reaches a stationary value when a certain 
characteristic value of Zc is reached. 

It is easy to estimate the characteristic value of Zc from energy considerations. The 
energy density of the first generation of particles in the magnetosphere of a magnetar cannot 
exceed the energy density of the primary particles accelerated in the inner gap. The energy 
density normalized to the magnetic field strength is 



3^3/6 /2 
32^ U 



a—R 

-Do 



1/3D2/3 (s) 4/3 



4/3- 



7min 
/max , 



Z 

Rs 



2/3 



(60) 



Given Eq. (60), the condition E[^'^^ < Nq^q/ {2Bq) transforms to 



7max 
7iiiin 



2/3 



< 



2a 



7n 



_P_ 

Ro 



100. 



(61) 



For a typical primary particle energy 70 ~ 10^—10^ the maximum energy of the 
first-generation particles is 7max — 10^—10^, Therefore, the generation of particles termi- 
nates at heights Zc — Rs — lORs- 
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An upper limit for the plasma multiplicity A = N^/Nqj in the magnetosphere of a 
magnetar can be easily obtained from Eq. (61): 

^max 

However, it should be noted that the electron-positron pair generation will be effective 
only if the inclination of the photon wave vector to the magnetic field direction is small 
compared to the maximum possible angle Xmax (see (18)). This is equivalent to the fact 
that the photon mean free path is much smaller than the distance from the dipole center 
to the point of photon emission, i.e., / ^ zq. If, however, / becomes equal to zq, then 
according to the factors mentioned above, the generation efficiency of a secondary plasma 
will be low. It is easy to find a criterion under the satisfaction of which the first generation 
of particles will be produced. For this purpose, we must equate the maximum possible mean 
free path If = p/7max to the distance from the stellar center zq = Rg. In this case, we must 
take into account the fact that the primary particle energy 70 appearing in the expression 

(c) 

for 7max is determined by the electric potential in the gap The potential distribution near 
the stellar surface can be obtained by solving the Poisson equation 

r dr dr dz'^ ) ^ ' 

where pe is the electric charge density and pcj = — r2B/(27rc) is the Goldreich-Julian charge 
density. The boundary conditions are a zero potential at the stellar surface z = Rg and at 
the surface defined by the last closed field lines, r = Rp, where Rp = R^J"^ / R^^"^ — 100 m is 
the polar cap radius. The solution is given in the form of a series in Bessel functions (Beskin 
et al. 1993). If the distance from the stellar surface {z — R^) > Rp/fii, where pi ^ 2.4 is 
the first root of the zeroth-order Bessel function, then the potential ceases to change along z 
and the expression for the radial potential distribution in dimensionless units is 

Bo R^. 



Rp 



Using the expression for 7max and given that 70 = |^'(r)|, it is easy to obtain the radial 
distribution of the mean free path for the photons emitted by primary particles near the 
stellar surface: 



8 R,.\ RsRu 



SBoRlJ (|cose|(l-zo))^^i 



1 - 



Rp 



(62) 



We see from Eq. (62) that the mean free path I reaches a maximum at r = Rp/2. 
Equating I to the stellar radius Rg, we can obtain a lower limit for the surface magnetic field 
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strength of a magnetar at which the generation of secondary particles is still possible: 

o p7/3 

5o>3^0, i?o>l, (63) 

28/3 



3| cos6'|(l — io) 

If the longitudinal current and the angle between the rotation axis and the magnetic dipole 
axis are small, then the coefficient (p ~ 2.1. Once the dimensions have been restored, we 
obtain 

/ p \ 7/3 / p \ -3 -^gi2 



DISCUSSION 

We may conclude that for each specific rotation period P, there exists some minimum 
surface magnetic field strength at which an effective generation of a secondary plasma begins. 
Qualitatively, this can be understood in the following way. The distribution function of 
electrons and positrons in the magnetosphere of a magnetar is limited below by 7min and 
above by 7max- As the radius of curvature of the field line increases, the minimum particle 
Lorentz factor increases, while the maximum one decreases. There exists some critical radius 
of curvature at which these become equal. Since the radius of curvature increases with 
neutron star rotation period, high primary particle energies proportional to the surface 
magnetic field strength are required for stars with long rotation periods. If we equate 'y^nm 
and 7max, while using a simple estimate of 70 — BqRI / {2Rl) for the primary particle energy, 
and take the radius of curvature at the point of intersection of the last closed field line with 
the stellar surface, then we will exactly obtain estimate (63) for the field strength but with 
the coefficient = 1. 

Equation (64) allows us to construct the line in the P—P diagram that separates the 
neutron stars with an active generation of electron-positron pairs in their magnetospheres. 
If we use the surface magnetic field Bq ~ 2(PP_i5)^/^ 10^^ G estimated from the formula 
for the magnetodipole losses and assume the longitudinal current and the angle between 
the dipole axis and the stellar rotation axis to be small, jy <^ jcj = cpcj and cos 6* ~ 1, 
respectively, then we can easily find that 

logP_i5 = y log P- 0.54, (65) 

where P is the neutron star rotation period in seconds and P_i5 is the period derivative 
in 10~^^ s s~^ 
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The figure presents the P—P diagram in which line (65) was drawn. When constructing 
the diagram, we used data from the ATNF catalog (Manchester et al. 2005). We see that 
all magnetars are located to the left of this line and fairly close to it. This suggests that 
secondary particles can be generated in principle in their magnetospheres. Therefore, we 
have no reason to believe that the generation of radio emission is suppressed in magnetars. 
Indeed, we showed here that the photon splitting, which was previously considered to be a 
pair generation suppression factor, is not such a factor. Although the splitting undoubtedly 
takes place in the superstrong magnetic fields of magnetars, it leads only to a strong po- 
larization of the curvature photons, but not to a decrease in their number or energy. The 
corresponding change in the maximum characteristic energy of the secondary particles can- 
not be a suppression factor either, because this energy for magnetars is even higher than that 
for radio pulsars. The ratios of the first-generation particle densities to the magnetic field 
strength in the cases of weak and strong fields are equal in order of magnitude. Therefore, 
the first-generation plasma density in the magnetosphere of a magnetar is even higher than 
that in themagnetosphere of a pulsar with a weak magnetic field. 

The fact that no secondary generation of particles is produced in a strong magnetic field 
can serve as a factor of suppression but incomplete one. Since the first-generation particles 
are produced straight at the zeroth Landau level, there is no emission of the synchropho- 
tons that produce the second-generation particles. The contribution from synchrophotons is 
known to dominate in the magnetospheres of radio pulsars. Numerical simulations of plasma 
generation cascades (Daugherty and Harding 1982; Hibschman and Arons 2001) suggest that 
the total secondary particle distribution function has the form of a power law -^n 
contrast, only the first generation of secondary particles with a distribution function propor- 
tional to 7~^/^ exists for magnetars. It thus follows that, in general, the plasma generation 
multiplicity A = N/Nqj in the magnetosphere of a magnetar is smaller than that in the 
magnetosphere of an ordinary radio pulsar. 

The already mentioned closeness of the minimum and maximum Lorentz factors for the 
first-generation particles, which can lead to a sharp decrease in the electron-positron plasma 
density (see Eq. (57)), can serve as another plasma generation suppression factor. Indeed, 
all magnetars are fairly close to line (65). In addition, the plasma generation efficiency 
depends on the longitudinal current flowing in the magnetosphere and on the angle between 
the axes. Obviously, this dependence is particularly sensitive to a change in parameters 
near the threshold line. The threshold magnetic field strength (64) also increases with 
current jy and angle 6. Therefore, adjacent magnetars in the P—P diagram may or may 
not exhibit an activity in the radio frequency band. For example, according to Malofeev 
et al. (2004), weak radio pulses are observed from AXP IE 2259+586 at a low frequency of 
111 MHz. According to Camilo et al. (2006), an intense radio emission at high frequency 
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0.7-42 GHz with a linear polarization close to 100% (Kramer et al. 2007; Camilo et al. 2007) 
emerged from AXP XTE 1810-197 after its X-ray outburst in 2003. The differences in the 
frequency and intensity of radio emission from these two magnetars could be explained by a 
difference in plasma density and multiplicity, since the first magnetar is considerably closer 
to the threshold line than the second one. However, no radio emission from other magnetars 
has been observed so far, which may be indicative of the influence of some other factors. 
Nonetheless, undoubtedly, one might expect the possible presence of weak radio pulses from 
them. 

For magnetars close to the threshold line (65) in the P—P diagram, the plasma density 
is low compared to that in the magnetospheres of ordinary pulsars. If, as is customary, the 
radio frequencies u are defined by the relation uj"^ ~ Up/'j^, where = AnNe'^/m is the 
square of the plasma frequency, then this leads to a shift in the observed radio frequencies 
toward the lower values, which is observed (Malofeev et al. 2005). 

We see that an effective electron-positron plasma generation is possible in the magne- 
tosphere of a magnetar. However, the absence of radio emission from most magnetars may 
stem from the fact that the presence of primary photons with a sufficiently high energy that 
must fall into the polar region is required to trigger a cascade pair generation. But its area 
is inversely proportional to the light-cylinder radius and, hence, is directly proportional to 
the rotation rate. At the same photon flux density from the external background, the prob- 
ability of triggering the plasma generation for magnetars with periods P ~ 10 s is a factor 
of 10—100 lower than that for pulsars with typical periods P ~ 0.3 s. We probably observed 
such a triggering in AXP XTE 1810-197. 
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logP 

Fig. 1. — Threshold hne defining the boundary of effective pair generation in the magneto- 
sphere of a magnetar in the P — P diagram. The crosses, squares, and dots mark SGRs, 
AXPs, and the remaining pulsars, respectively. For comparison, the death line for ordinary 
pulsars with a slope of 11/4 is also drawn. 



